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Abstract. We renormalize QCD at three loops in the modified regularization invariant, RI', 
scheme in arbitrary covariant gauge and deduce that the four loop /3-function is equivalent 
to the MS result. The anomalous dimensions of the scalar, vector and tensor currents are 
then determined in the RI' scheme at three loops by considering the insertion of the operator 
in a quark two-point function. The expression for the scalar current agrees with the quark 
mass anomalous dimension and we deduce an expression for the four loop RI' mass anomalous 
dimension in arbitrary covariant gauge and for any Lie group. 
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1 Introduction. 



The most widely used renormalization prescription in perturbative quantum field theory is the 
minimal subtraction scheme where only the infinities with respect to the regularization are 
subtracted from the divergent part of the Green's function to determine the renormalization 
group functions, Jj. In practice, though, it is more appropriate to use the modified minimal 
subtraction scheme, MS, since the convergence properties of perturbative series in this scheme are 
improved by additionally absorbing a finite part, ln(47re~'^), into the renormalization constants, 
[2], where 7 is the Euler-Mascheroni constant. The advantage of using the MS scheme, which is 
a mass independent scheme, rests in some elegant properties. For example, the MS /9-function 
and anomalous dimension of the quark mass in QCD, are both independent of the covariant 
gauge fixing parameter, ^Ej- Moreover, performing computations using MS and dimensional 
regularization, where the spacetime dimension becomes d = 4 — 2e and e is the regularizing 
parameter, one can carry out multiloop calculations to very high order. Indeed in QCD various 
four loop renormalization group functions are available, [^[3115], which represent the current 
state of computation. Whilst the MS scheme enjoys these elegant features and has become 
the standard reference scheme, it has the limitation that it is not a physical renormalization 
scheme. Examples of schemes which are founded in a more physical origin include, for example, 
the MOM and MOM schemes, |71 |H]. Though one disadvantage of using physical schemes is 
that fewer results currently exist to the same multiloop precision as in MS. However, it is 
well known that physical quantities in one scheme can be simply related to the same quantity in 
other schemes by a conversion function, [Jj. Whilst one ordinarily uses dimensionally regularized 
perturbation theory to compute physical quantities one can also determine such information by 
using a lattice regularization. The advantage of this approach is that one in principle includes 
all non-perturbative contributions in a calculation which need to be converted from the lattice 
scheme to MS. For a recent review and applications in determining matrix elements in deep 
inelastic scattering see, for example, (2j where lattice results were matched to MS results. The 
scheme used is similar to a modified version of the regularization invariant, RI, scheme known as 
the RI' scheme, jTH]. Therefore, in order to improve lattice estimates one requires the conversion 
of various renormalization group functions from MS to RI'. Recently, this problem has been 
addressed in the context of quark masses where the conversion functions were produced for 
all covariant gauges for the quark mass anomalous dimension at four loops, jlll I12j . Indeed 
the field anomalous dimensions were also deduced to the same order for the SU{Nc) colour 
group, Though in practice for the lattice application one only considers one particular 

gauge which is the Landau gauge. This work of J2] extended the three loop calculation of 
jllj . One practical feature of these computations was that one only needed to consider ordinary 
perturbation theory in the massless limit which effectively meant that the conversion functions 
could be deduced using standard multiloop perturbative tools for massless field theories. 

Whilst these papers dealt with the problem of quark masses deduced from the lattice there 
are other problems where the conversion functions are required. For instance, there is interest 
in deducing low moments of the structure functions measured in deep inelastic scattering from 
the lattice, [HI El- To improve estimates the conversion factors from the MS scheme to the 
RI' scheme are required. Therefore, the purpose of this article is twofold. First, given that 
the RI' scheme is important for relating Landau gauge lattice results to the MS scheme we 
will renormalize QCD at three loops in the RI' scheme though in a general covariant gauge. 
Whilst we are ultimately interested in quark currents it is not inconceivable that the anomalous 
dimensions of operators with gluonic fields will at some time be measured on the lattice and 
therefore the anomalous dimensions of the gluon (and ghost) fields will need to be determined at 
the same level as the quarks. Equipped with the fully renormalized QCD Lagrangian in the RI' 
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scheme we will then extend the approach of ^21 to deduce RI' information but for the anomalous 
dimension of a particular quark composite operator which corresponds to the tensor current in 
QCD. It is of interest since it represents the lowest moment of the transversity operator in deep 
inelastic scattering, ^1]. Given the recent resurgence of experimental and theoretical interest 
in transversity, (see, for example, ^| El ^1 El El)) the longer term aim is to provide a 
more accurate numerical estimate for the associated matrix element prior to experimental data 
being accumulated at RHIC. One motivation in this approach is to develop the calculational 
formalism to determine the conversion function for an operator which is a simple extension of 
|12j . In |12j the quark mass conversion functions were determined by considering the corrections 
to the massive quark two-point function which is an avenue not immediately available for a 
composite operator. Nevertheless we will bridge this gap by reconstructing the result of ^21 at 
three loops by first considering the problem of the mass renormalization as the renormalization 
of the associated composite operator, ipip, as inserted in a two-point Green's function prior to 
replacing it by the operator of main interest which is the tensor current. Moreover, given this 
way of computing we need only use the massless version of QCD. 

The paper is organised as follows. In section two we discuss the three loop renormalization 
of QCD in the RI' scheme and provide the renormalization group functions. These results are 
used in section three to extract the anomalous dimension of the quark mass operator in QCD 
in RI' which agrees with the earlier Landau gauge result of |12j . Having provided the formalism 
for treating an operator, we extend that calculation to the tensor current case in section four. 
Finally, our conclusions are given in section five. 



2 RF scheme at three loops. 

We begin by explicitly renormalizing QCD in the RI' scheme at three loops using dimensional 
regularization. The bare QCD Lagrangian, with the gauge fixed covariantly, is 

L = - Ig" G"'^'' - — - c''^^'D^c'' + 0^ Ipil,^^ (2.1) 

where A'^^ is the gluon field, G^j, = dfj.A'^ — d^A'^^^ — gf"''"'A^^A'^, c" and c"- are respectively the 
ghost and antighost fields and a is the covariant gauge fixing parameter. The indices range over 
1 < a < Na, 1 < I < Nf and 1 < i < Nj where Np and Na are the respective dimensions of 
the fundamental and adjoint representations of the colour group whose generators are and 
structure functions are /"'"^ whilst Nf is the number of quark flavours. The covariant derivatives 
are defined by 

D^^P = d^i, + igA^T'^i^ , D^Al = d^Al - gr'^A^Al 

D^c^ = d^c^ - gr^'Alc' . (2.2) 

In renormalizing the full Lagrangian in the RI' scheme, which we will define more precisely 
later, we must ensure that the scheme is consistent. We can achieve this, for example, by 
demonstrating that the renormalization constants for the gluon, ghost and quark fields correctly 
produce the same gauge independent coupling constant renormalization at three loops from 
different Green's functions, thereby ensuring that the Slavnov- Taylor identities are respected. 
If we regard all the quantities in the QCD Lagrangian, (|2.1() . as bare and denote them with the 
subscript o, we introduce the corresponding renormalized quantities by the usual definitions 

^7 = ^A^'' , = ^/Z'.c" , Vo = V^V' , go = f^'Zgg , ao = ^"^^a a (2.3) 
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where /x is the mass scale introduced to ensure the couphng constant is dimensionless in d- 
dimensions and d = 4 — 2e with e the regularizing parameter. 

To determine the RI' scheme values of the renormalization constants we first consider the 
gluon, quark and ghost two-point functions. In IIH I12j the RI' scheme definition of the 
quark wave function renormalization is given by the Minkowski space condition 



lim 



(2.4) 



where S^(p) is the bare (massless) quark two-point function and p is the external quark mo- 
mentum. As we are considering massless quarks S^(p) will be proportional to and involve 
poles in e at each order in the strong coupling constant. To contrast with the MS scheme, the 
RI' scheme definition of the quark wave function is such that one absorbs the complete finite 
part of the Green's function with respect to e into the renormalization constant. In other words 
only the 0(1) piece is removed and the 0(e) part is ignored. In the MS scheme only the poles 
in e are removed as well as the finite parts involving powers of ln(47re~'*'). For completeness we 
note that in the RI scheme one absorbs the full finite part in the same way as in the RI' scheme, 
(|2.4|1 . but for a different part of the Green's function which is. 



lim 



1 

4d 



Z. 



RI 



d 



1 



(2.5) 



Due to the presence of the derivative this scheme is much more difficult to implement on the 
lattice compared with H2.4() which is why we are concentrating on RI'. For the remaining wave 
function renormalization constants we define their RI' values in a similar way to (|2.4I) . For the 
ghost fields is determined by 



lim 



z 



RI' 



(2.6) 



The definition of the gluon field requires more care due to it having transverse and longitudi- 
nal components. Rendering the former finite determines Za whilst the finiteness of the latter 
component fixes the gauge parameter renormalization constant. In the spirit of the quark and 
ghost RI' scheme renormalization constants we now define Z^^ and Z^^ through the following 
conditions. Writing the gluon polarization tensor, Il^^{p), as 



n^.(p) 



Urip) 



p2 



+ ^l{p)j^ 



(2.7) 



we define the associated RI' scheme renormalization constants as 



and 



lim 



lim 



zfUrip) 



ZfULip) 



(2.8) 
(2.9) 



Having established the scheme definition of the field renormalization constants we have com- 
puted them explicitly for arbitrary a for 1)2. 1(1 . We have used the Mincer package, [20], written 
in the symbolic manipulation language Form, |211l22j . where the Feynman diagrams are gener- 
ated with Qgraf, |23]- The basic nature of the renormalization conditions are straightforward 
to implement in FORM using the approach of j24j . Briefly, one computes the appropriate Green's 
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functions in terms of bare parameters and couplings and then introduces the counterterms by 
rescahng these variables into the renormalized ones using the definition of renormalization con- 
stants, (|2.H|) . Also, for Green's functions which ordinarily involve a bare parameter in the tree 
part, one divides the Green's function by the bare parameter first before rescaling, [2^. There- 
fore, we quote the results of our calculations. We find 
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where Tr (t'^T^) = ^p^°^ TT" = Cpl, f'""^f^'"^ = CaS"^, C{n) is the Riemann zeta function, 

a = 5^/(167r^) and g is the couphng constant appearing in the covariant derivative. At this 
point it is worth commenting on the status of the variables a and a. As we have computed 
the renormahzation constants for the RI' scheme they correspond to the RI' scheme couphng 
constant and covariant gauge parameter. When it is necessary we will include a label on the 
variables to distinguish which scheme they are defined in. As with other schemes they can be 
related to the corresponding MS variables which we will discuss later. Throughout this article 
when wc quote any renormahzation constant the scheme will be denoted on the renormahzation 
constant itself and it will be understood that the variables will be in that scheme as well. For 
completeness we note, 
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where the variables a and a correspond to those of the MS scheme. There are several checks on 
these results. First, we have verified that the correct three loop MS renormalization constants, 
[251 1261 1271 1281 129j , emerge with the programmes we have written prior to extracting the results 
for the RI' scheme. Second, the gauge parameter, a, does not get renormalized as in the MS 
scheme so that gauge invariance is not destroyed. Next, our three loop result for agrees 
with the four loop result of ^2] in the Landau gauge. They do not agree for non-zero a. The 
reason for this is that in the RI' scheme used in I12J the values taken for Za and Za were different 
from those given above. More specifically the MS expressions were used which only differ in the 
case a 7^ 0. However, using the scheme adopted in we have reproduced the expressions given 
there for all a which again confirms the correctness of our programming. However, as we believe 
the full RI' scheme we have introduced here is more natural we will always present subsequent 
results with reference to the above renormalization constants. Indeed it is not inconceivable 
that at some point one would require the renormalization of an operator with gluon content in 
RI' and therefore one would require our Z^ . Nevertheless all results should be consistent with 
[12] in the Landau gauge, q = 0. 



Whilst the results for the field renormalizations are deduced from the two-point functions 
we have yet to establish their consistency. This can be verified by renormalizing several three- 
point functions to obtain the same coupling constant renormalizations. We have examined both 
the quark gluon and ghost gluon vertices. However, with the wave function renormalization 
containing a finite part at one loop the method of renormalizing the vertices in the RI' scheme 
cannot be the same as for the two-point functions. In other words one cannot absorb a finite 
part from the three-point functions into the definition of the coupling constant renormalization. 
For instance, at one loop the quark gluon vertex would give an Nf dependent finite part to that 
coupling constant renormalization constant. However, this dependence cannot be matched 
in the gluon ghost vertex which is independent at one loop. Although it is possible to follow 
this route and accommodate the problem of the finite parts in the three-point vertices it will 
lead to an MOM or MOM renormalization scheme which we are not considering here, [ZllHl- 
Therefore, to define a three loop renormalization constant for Zg in the RI' scheme which is 
consistent with the Slavnov- Taylor identities we define the Zg renormalization for each vertex 
in an MS way. In other words we define Zg by only absorbing the infinities without removing 
any finite parts, aside from powers of ln(47re~'''). With this we can consistently deduce the same 
coupling constant renormalization from both three-point functions. More concretely we have 
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and 
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where Gf^^^*^ are treated as finite with respect to e and are not unity and we have omitted the 

overall structure function from the ghost gluon Green's function. Though the tree part of each 
(^finite ^^j^j^y 

since we have first divided Gi{p) by the bare coupling constant before rescaling 
as in the method of To evaluate Z^^ explicitly we have again used Mincer, ^20,. However, 
to do this correctly the external momentum of one leg must be nullified as Mincer can only 
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be applied to two-point functions. To avoid potential spurious infrared infinities arising in this 
case we have chosen to nullify an external quark or ghost leg respectively leaving p as the overall 
external momentum. Consequently we find from both three-point functions that 
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which is the same as in the MS scheme, |291 128| . though it could only have differed in the three 
loop term which is scheme dependent. Thus to this order the RI' and MS /3-functions coincide. 
In our conventions we have, for all a. 
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where we have explicitly indicated the scheme of the coupling constant as a subscript. Given 
that the three loop term of the /3-function can be different we have constructed the relation 
between the coupling constants of the RI' and MS schemes explicitly. To do this we parallel 
the approach of (30j where the same procedure was followed to establish the relation between 
various MOM coupling constants and the MS coupling. Since we have renormalized two different 
three-point functions we have to check that both give equivalent results and are consistent with 
the Slavnov- Taylor identities. We define 
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where the Green's functions on the right hand side are the finite expressions. Moreover, since we 
are concerned with the finite part of the Green's function after renormalization G^}\ corresponds 



to a particular Lorentz projection of the full Green's function G^^^, 
momentum routing we are considering, we have defined, j30j . 



PU]. In particular, for the 
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where 



(2.22) 
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For the ghost gluon vertex we do not need to take into account several projections since there 
is only one Lorentz structure for that vertex with one external momentum nullified. In addition 
to computing as a function of for both vertices we will also need the relation of the 

covariant gauge parameter in one scheme with that in the other since, for instance, G^g^^"^^^® 
depends on a-^ and aj^j'. This is achieved by, |3nj . 

^RI' 

"RT' = ^^"aIq • (2-23) 
7MS MS ! 

^A 

By solving (|2.19j) . (|2.20j) and (|2.23|) iteratively we have determined the relationships between 
the coupling constants and covariant gauge parameters in both schemes at three loops. Prom 
both the vertices we find 
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- 884736C(3) - 1343872) CATjrNf 
+ ( - 3068928C(3) + 4976640C(5) - 988416) Cj^TpNf 

+ (2101248C(3) - 2842368) CpTj^N?^ 



31104 



"MS + O ay . (2.25) 



The former expression is consistent with the three loop RI' /3-function being equivalent to the 
MS one for arbitrary covariant gauge. Moreover, since the three loop term of the transformation 
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is also absent this implies that the four loop RI' /3-function is also equivalent to the four loop MS 
/3-function in all gauges. The non-trivial relation between the gauge parameters will be crucial 
in carrying out checks on the renormalization group functions. 

We have calculated the renormalization group functions for the various wave function rcnor- 
malizations directly from the renormalization constants themselves. In particular we used 



ifia) 



7?' (a) 
ifia) 



+ «7a (a)- 



(5{a 
/3(a) 



da 

da 



da 



a- 



dlnZ^ 
da 



+ a7^^ (a)- 



da 
dlnZf-^' 
da 



+ «75(«) 



da 
dlnZ^'^' 
da 



though Z^^ = 1 at three loops implies that 



(2.26) 



(2.27) 



to the same order which corresponds to the gluon propagator being transverse. Hence, from the 
renormalization constants we have computed we find, in four dimensions, that 



jf{a) = [STpNf - {IS - 3a)CA] 



6 



(27a^ - DOa^ - 426a + 372?) Cl + {l2a^ + 240a - 3616) CaTfN^ 



hVimrlNf - 15552ClTFNf + (331776C(3) - 487296) CpT^N^ 

- (486a^ + 3078a^ + 10260a3 - 1458C(3)a2 - 25965a2 + 86184C(3)a 

- 173406a - 175446C(3) + 2127823) Cl - (648a^ + 216a^ + 47808a2 
+ 10368C(3)a + 126480a - 254016C(3) - 2501184) ClTpNf 

- (7776a2 - 62208C(3)a + 71280a + 725760C(3) - 1131408) CaCfTpN} 



216 



+ 



+ (^11520a^ + 19200a - 165888C(3) - 751680j CaTfN/ 



7776 



+ O(a^) 



(2.28) 



7f' (a) 



aCpa + 
+ 



(da^ + 45a^ + 223a + 225) Ca - 54Cf - (80a + 72) TpNf] 

(l62a^ + 1377a'^ + 7578a^ + 1134C(3)a2 + 22608a2 

- 23004C(3)a + 113080a - 39690C(3) + 179811) Cl 

- (648a^ + 1944a2 - 15552C(3) + 52272) CaCf 

- (l440a^ + 7200a2 + 5184C(3)a + 63616a - 10368C(3) + 110016) CaTeN} 
+ (20736C(3)a - 23760a + 6048) CpTrNf 



+ (6400a + 14976) T^Nf + 1944C|. 



1296 



+ 0(a^ 



(2.29) 



12 



and 



[a - 3J — — + 



+ 18q^ 



813] Ca - (80q - 312) TpNf 



144 



+ 



(l62a^ + 891q^ + 972C(3)q^ + 1503a^ + 4050C(3)q2 - 20700^ 
- 15876C(3)a + 46363a + 34182C(3) - 471909) Cl 
- (l440a^ + 2880a2 + 7776C(3)a + 39208a - 33696C(3) - 32268o) CaTfNj 
+ (20736C(3)a - 23760a - 62208C(3) + 79056) CpTpN} 

■ (2.30) 



+ (6400a - 48000) T^N^ 



5184 



+ O(a^) 



where we use the same convention for the renormahzation group functions as for the renormal- 
ization constants in that the variables a and a are in the scheme indicated on the renormahzation 
group function itself. The expression for (a) agrees with the three loop expression for the 
colour group SU{Nc) in the Landau gauge given in Moreover, a final check on our calcu- 
lation resides in the fact that in constructing these RI' scheme renormahzation group functions 
the correct double and triple poles in e in the renormahzation constants have been determined 
in the computation. If they were not correct then finite renormahzation group functions would 
not have emerged. 

In ^2] the quark anomalous dimension was computed explicitly by first determining the 
appropriate function which converts the MS result to the RI' expression following a standard 
procedure which is discussed in, for example, As a final check on our wave function renor- 
mahzation group functions in the RI' scheme we have also computed them from the conversion 
functions which are defined as 



^RI' 

CA(a,a) = 



Cc(a, a) 



z. 



RI' 



1= , C^(a,a) = -= . 



ZMS 



Z^ 



(2.31) 



It is important to appreciate how these functions are explicitly constructed. They are functions 
of the two parameters a and a in the same scheme. However, the RI' scheme renormahzation 
constants depend on Oj^j/ and a-^ji which therefore must be converted to their MS counterparts. 
In the following expressions for the conversion functions, and those we give later, we have chosen 
to express them in terms of the MS variables and omitted the corresponding subscript. We found 



CA{a,a) 



1 + 



(da^ + 18a + 97) Ca - SOTpNf 



a 
36 



+ 



(sWa^ + 2430a2 + 5184C(3)a + 2817a - 7776C(3) + 83105) Cj 



- (2880a + 20736C(3) + 69272) CaTfNj + (41472C(3) 
+ 12800T|A^/ 



47520) CpTpNf 



2592 



+ 



17010C(5)a^ - 12636C(3)a'^ + 64638a'' - 51516C(3)a^ + 19440C(5)a^ 
+ 322947a3 + 203148C(3)a2 - 8748C(4)a2 - 374220C(5)a2 + 1094553a2 
+ 4636764C(3)a - 34992C(4)a - 1341360C(5)a + 1457685a 
- 7531056C(3) - 26244C(4) - 3414150C(5) + 44961125) Cl 
+ (l5552C(3)a2 - 303912a2 - 3670272C(3)a - 890064a - 7293888C(3) 

+ 839808C(4) + 4976640C(5) - 49928712) ClTpNf 
+ (746496C(3)a - 855360a + 23514624C(3) - 1119744C(4) 



13 



+ 
+ 
+ 



+ 7464960C(5) - 37099872) CaCfTfNj 
(663552C(3)a + 64512a + 5971968C(3) + 13873536) CaTInJ 
(9206784C(3) - 14929920C(5) + 2965248) C^TpNf 



(16130304 - 12939264C(3)) CpT^Nf 



1024000r|Af/ 



93312 



O(a^) 



(2.32) 



Cc(a, a) 



1 + CAa 



+ 
+ 



72a" 



CAa^ 
192 



( 29241a 



36C(3)a2 + 72C(3)a - 21a - 180C(3) + 1943j Ca - imTpNf 
^ - 11178C(3)a3 - 4860C(5)a3 - 56862C(3)a2 + 1458C(4)a2 

+ 34020C(5)a + 1027890^ + 254826C(3)a + 5832C(4)a - 63180C(5)q 
- 3510a - 728082C(3) + 4374C(4) - 63180C(5) + 4329412) C\ 

+ (42984Q - 67392C(3)a - 100224C(3) - 139968C(4) - 2650192) CaTfNj 

+ (684288C(3) + 186624C(4) - 1165104) CpTpN} 



+ (27648C(3) + 330304) T^Nf 



CAa'' 
31104 



+ 0{a^) 



C^{a,a) 



+ 



+ 



aCpa 

Ba^ + 5) Cf - (90^ - 24C(3)a + 52q - 24C(3) + 82) Ca + 28TFNf 

(l728C(3)a^ - 11880a^ + 25272C(3)a2 - 972C(4)a2 - 6480C(5)a2 

- 63747a2 + 181440C(3)a - 1944C(4)a - 12960C(5)a - 358191a 
+ 678024C(3) + 22356C(4) - 213840C(5) - 1274056) 
+ (123120^ - 5184C(3)a^ - 31104C(3)a2 + 596160^ + 57024C(3)a 
- 103680C(5)a + 85536a - 228096C(3) - 31104C(4) 
+ 103680C(5) + 215352) CaCf 
+ (3456C(3)a^ - 5184a^ - 11016a + 31536) C| 
+ (124056a - 41472C(3)a - 89856C(3) + 760768) CATpNf 
+ (68256 - 82944C(3) - 28512a) CfTeNj 
Cpa^ 



(2.33) 



Cpa^ 



100480r|A^/ 



5184 



+ 0(a^ 



(2.34) 



With these the RI' scheme renormahzation group functions can be determined from 



^RI' 



{a 



RI' 



MS 



(""Ms) + ^ ("ms) ("^Ms ' "ms, 

IVIS 



Q 

+ "ms'''" i°Ms j f)n, v"ms ' "Ms j 

aa^g 



(2.35) 



where i = A, c or tp and we have included the scheme dependence of the variables explicitly 
though to the order we are working to the /3-function is the same in both schemes. For each of the 
three cases we have computed the right hand side of (|2.35|) in terms of the MS variables and then 
converted to their RI' counterparts before verifying that the same previous expressions correctly 
emerge in terms of the RI' scheme variables. This completes the full three loop renormahzation 
of the QCD Lagrangian in the RI' scheme. 
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3 Quark mass anomalous dimension in the RI' scheme. 



We now turn to the problem of deducing similar renormalization constants for the composite 
quark currents of the form 0_4 = t/^Ail^ where ^ = 1, 7^ or a^'^ with a'^'^ = ^[7^,7'^] where the 
latter corresponds to the tensor current. As the lattice data is available for the insertion of the 
operator at zero momentum to deduce the conversion functions we will insert each of at zero 
momentum into a quark two-point function and examine the divergence structure of 



GoAp) = i^ip) [Mm) H-p)) = {^{p) Oa{^) i>{-p)) 



(3.1) 



In order to demonstrate the validity of this approach we must first reconstruct the anomalous 
dimension for the mass, ^^1) which corresponds to the operator ^ = 1. Therefore, defining the 
renormalization constant Z^^ of the operator in the usual way by, 



the RI' scheme value is given by the condition 



(3.2) 



lim 



(3.3) 



where the wave function renormalization constants arise from the external fields. In the MS 
scheme this results in a gauge independent renormalization constant to all orders, However, 
in RI' this is not the case since 



Z 



RI' 



1 



+ 



+ 



+ 



+ 



+ 4 + a Cpa + 



11 



-Cp + —CfCa - 2TFNfCF 



1 



5 97 
-TpNfCp — —CfCa + 



45 



83 20a 



19 



TfNjCf + 



4 

18C(3) 



+ 3a ) 

1285 



223a 5a' 



+ — - 12C(3) +4a + a^ 



16, 



24 

,2 



a 



36 



121 



a 
T 



33 



CfCa 



—TfNjCfCa + QTpNfC'p - -T^N/Cp - —CpCi - —C'pCA - -C'p 



+ 

+ 
/556 



^^T^NfCp - ^TpCpCA + ( 2a 



27 
49 
12 



11a 



27 

CI'Ca 



TpNjCl + ^^CfC\ 



63 9a 

T + T 



16C(3) ) TpN^CfCa - + 16C(3) + ^ 



o 140 o n 

TpNfCl + —T^NfCF 



11413 2 ^ / 4889 



r.^/oN 80a 15a2 'io?\ ^1 ^ 
54C(3) + — + ^ + — CICa 



+ 



616 



C(3) - 4aC(3) 



4 

95387 
243 



cl\- 



+ [^C(3)+8aC(3)+24C(4) 



24C(4) 

1109 
~9" 



3976a 50a' 



lOa'^ 



81 



9 



TpNfC fCa 



+ 



115a 40a2 



18 



+ 



TpNfCl 



15 



+ 



/ 7514 , 400a , 32 
V 243 



/ 3360023 28351a 



157a2 



\^ 3888 



324 



42 la'^ ITa'^ a^ 
+ + —r^ + 



+ 



72 
/962 



16 



31193 
72 



C(3) 



+ 43a + 3a^ C(3) - 20C(5) 



^C(3) + ^C(3) + 60C(5) j CpCl 
18781 6089a 



72 



72 



170a' 



4a^ 



a 



, C%C A + I 
2 / ^ V 12 



/3227 31a 



+ 



58 - 6a + 6a' C(3) - 120C(5) C% 



+ 4a' + a-^ 



(3.4) 



which agrees with the quark mass renormahzation constant deduced at four loops in 
in the Landau gauge in our conventions. Thus we have demonstrated the equivalence of our 
operator method with the massive propagator approach of jl2| . Further, we have checked 
that the correct three loop MS quark mass anomalous dimension, [321 133j . emerges from our 
programmes. Therefore, from we can deduce the corresponding renormahzation group 

function. With 

ainz^f 

- f3{a) 



7g'(a) 



da 



ai^' (a) 



da 



(3.5) 



we find 



Cpa^ 



+ 



SCpa - [(185 + 9a + 3a2)Cyi + 9Cf - 52TFNf] 
(l08a^ + 324a2 - 1944 - 19008C(3)) CaCf 
- (l 17428 + 5634a + 19050^ + 405a^ + 54a'^ - 28512C(3)) Cl 
(480a2 + 2088a + 6296o) CaTfNj - 13932C| 

+ (16632 - 3456C(3)) CFTpiV/ - 6848r|iV/) ^ 



+ 



216 



+ 0(a4) (3.6) 



which agrees with J2] in the Landau gauge, aside from an overall factor stemming from our 
conventions which are the same as jMl ■ We have also derived the same expression by constructing 
the conversion function (a, a) with A = 1 where 



Co^{a, a) 



7RI' 



(3.7) 



Thus, using 



7o 



RI' 



Msj 



\ d 



"ms7^^ 



'MS 



Msj 



MS] Q^. 



MS 



(3.8) 



with the explicit expression 

C^^(a,a) = 1 - (a + 4)CFa + [(24a2 + 96a - 288C(3) + 57) 
+ 332rFiV/ - (l8a^ + 84a - 432C(3) + 1285) Ca 



Cfo^ 



24 



16 



+ 



+ 



(l5552a3 + 89424q2 - 23328a2C(3) + 573804a 
- 334368aC(3) - 2493504C(3) + 155520C(5) + 2028348) CaCf 

- (l3122a3 - 8748a2c(3) + 716850^ - 103032aC(3) 

+ 357777a - 3368844C(3) + 466560C(5) + 6720046) C\ 
113400a - 31104C(3)a 

- 532224C(3) + 186624C(4) + 3052384) CATpNf 
+ (62208aC(3) - 123120a - 331776C(3) 

- 186624C(4) + 958176) CfTfNj 

- (7776a2 - 46656a2c(3) + 311040^ + 46656aC(3) + 30132a 

- 451008C(3) - 933120C(5) + 2091096) C| 

Cpa^ 



(27648C(3) + 240448) T/A^. 



'/ 



7776 



+ 0(a*) 



(3.9) 



we find exact agreement. Moreover, this conversion function agrees with that given in JLSj when 
restricted to the Landau gauge. However, in checking the three loop expression in the RI' scheme 
only the contribution up to and including the two loop term is required. The three loop part is 
only relevant for the four loop anomalous dimension. Therefore, since the four loop MS quark 
mass anomalous dimension is available, 01^1) for an arbitrary colour group we can deduce that 
the four loop correction to 1)3. 6() is 



7g'(a) 



3Cfo - [(185 + 9a + 3a^)CA + 9Cf - 52Tp7Vf]— ^ 

6 

(l08a^ + 324a2 - 1944 - 19008C(3)) CaCf 
- (ll7428 + 5634a + 1905a2 + 405a^ + 54a'' - 28512C(3)) Cl 
+ (480a^ + 2088a + 6296o) CaTfNj - 13932C| 

+ (16632 - 3456C(3)) CfTfNj - 6848r|iV/)J 

1215a6 - 13608a^ - 90801a^ - 8262C(3)a^ - 368064a3 + 104004C(3)a2 

■ 1397826a2 + 940734C(3)a - 4554684a + 39004740C(3) 
1710720C(5) - 92569118) C^Cf 

(2916a^ + 28188a^ - 23328C(3)a^ + 136404a^ - 252720C(3)a2 

+ 3771360^ - 1717200C(3)a + 429300a - 22203072C(3) 

- 1710720C(5) + 10355148) C^Cl 
(l2960a^ + 103032a^ + 34992C(3)a2 

+ 677952a2 - 316224C(3)a + 3021840a - 14239152C(3) 

- 1244160C(5) + 73217928) ClCpTpA^/- 

- 3888a^ + 46656C(3)a^ - 11664a^ + 241056C(3)a2 - 5832a2 

- 1236384C(3)a - 103032a - 1601856C(3) 
+ 10264320C(5) - 33960384) CaC^^ 

+ ( - 25920a3 - 62208C(3)a2 + 18144a2 + 694656C(3)a + 230688a 
+ 1347840C(3) - 1244160C(5) + 20983248) CAClTFNf 



+ 



+ 



17 



+ (- 576000^ + 82944C(3)a - 449280a 

+ 580608C(3) - 16599552) CaCfT^N^ 

+ (- 62208a2 + 373248C(3)a + 31104a 

- 3856896C(3) + 9745920) C^TpNf 
+ (- 165888C(3)a + 41472a + 2571264C(3) 



6653952) C]pT^N? 



+ (2612736C(3) + 1225692) Cj, + 1025536CFr|7V/ 

Mbcd^abcd 

+ (248832 - 1866240C(3)) ^ ^ 



+ (3732480C(3) - 497664) Nj 



jiahcd^ahcd 
F F 



7776 



+ 0(a^ 



(3.10) 



^abcd ^abcd ^abcd ^abcd 

where ^ j^J" — and ^ j^J^ — are the quartic Casimirs associated with light-by-light topologies, 
[Uin], and the coupling constant and gauge parameter are in the RI' scheme. The four loop 



expression is in agreement with the Landau gauge expression of J2] for an SU{Nc) colour group. 

As an additional check on the renormalization of a composite operator in the RI' scheme we 
have also considered the case of „4 = 7^. One reason for examining this operator arises from the 
fact that as it is now a Lorentz vector the Green's function, G'^^^^(p), is not merely proportional 
to j^. Instead by Lorentz symmetry 



+ s 



(_2) 



ip) 



p2 



(3.11) 



where the amplitudes S^^,j^(p) depend on the coupling constant. They are determined by the 
relations 



tr 



4{d-l) 



{l^.G'^^,^{p) 



p 



The renormalization constant for the operator is determined from 



and renormalizing the Green's function we find 



Z 



MS 



1 + O(a^) 



(3.12) 



(3.13) 



(3.14) 



The non-renormalization of this current rests in the fact that it corresponds to a physical operator 
and therefore on general grounds its anomalous dimension vanishes. (See, for example, |31j.^ 
Furthermore, as the vector current has been inserted at zero momentum the 7^ component of 
its Green's function must obey the Slavnov- Taylor identity and be equivalent to the finite part 
of the quark two-point function after renormalization in the same scheme. In computing the 
quark wave function anomalous dimension in the previous section we have also determined the 
finite part in the MS scheme and it is reassuring to note that both it and 



(1) MS finite. ^ 



1 + aCpa 



+ 




3(l + a)C(3) CfCa 



18 



7 ^2 

-TfNjCf - -Cp 



a 



+ 



6 2 

52^, , , 11887 1723a 

+ 1 — C(3) + 8aC(3) 



81 TlAT/ + h 6C(3) - ^ - ^ ) TpNfCF 



81 

/ 159257 39799a 787a2 



\^ 648 ■ 576 
/ 3139 



24 



+ 35a + 



64 
39a2 



72 

55a^ 
24 



TpNfCA 



a 



+ y a3) 



3a^ 



+ 



3a _ 69 
~ 16 



C(4) 



/ 165 5a 5a^\ , 

+ 20(a - 1)C(5) + 6C(4) + (44 - 17a + a^) C(3) 

997 , 3a2 a3\ ^ ^ 

h 4a H CfCa 

24 2 8 / 



C^a^ + O(a^) (3.15) 



are in exact agreement which provides an additional check on our programming. Further, it is 
the finite part of the first term which determines the relation of the MS scheme to the RI' scheme. 
Such a feature of extra contributions will persist in the tensor current case but the vector operator 
is peculiar in the sequence of dimension three operators in that only its anomalous dimension 
and finite part are entwined with the Slavnov- Taylor identity. Repeating the same exercise for 
the RI' scheme by introducing the definition 



1 



(3.16) 



we find 



7RI' 



1 + O(a^) 



(3.17) 



This is consistent with the observation that if an anomalous dimension of a physical operator 
vanishes in one scheme it vanishes in any other scheme. Moreover, (|3.16|) is also consistent with 
the Slavnov- Taylor identity in the RI' scheme since not only is it finite prior to renormalization 
but it is also unity which agrees with the finite part of the quark two-point function consistent 
with the nature of this renormalization scheme. Finally, to assist with lattice matching we record 
that the finite parts of the second component of the vector current Green's function are 



^(2) MS finite, X 



2aCFa 



Cf 
Cf 



— + 7a + -a^ J Ca - 4rpA^j + {2a'^ - 3) Cf 



28 \ ^ . 208^2 ^^2 . 17 



— - llaj TfNjCf + -^T^Nf 



f»]C'F 



+ 



16C(3) 8C(3)a 



1528 175 \^ 



19 



/ 19979 4393 295 o 27 , 

+ 1 a H a H 

V 72 48 16 8 



/245 59 9 o\ 

/ / o\ 242 

+ ((24 - 6a - Ga^j C(3) - — + 33a 



+ 17a2 + ^a-" ) CfCa 



(3.18) 



and 

^(2) RI' finite, X 



— 2aCi.-a 

- Cf 



- Cf 



25 223 5 ^ 1 o\ ^ 
2 18 2 2 y 
40 \ 

4 + y a j TpiV/ - 3Cf 
28 110 

y ~ ~3~ 



a + 32C(3)a^ TfNjCf + 
+ (^16C(3) - 8C(3)a - 



/ 208 800 ^ 2 ;v2 



1528 7952 



a 



+ 



81 

100 2 20 

/ 19979 14135 314 2 421 , 17 4 1 . 

+ -sT^ + ^« + + + 4" 

- + y « - la^) C(3)) Cl + 3C| 



242 

+ ( 24C(3) - _ - 3a2 - a^ ) CfCa 



+ O(a^) (3.19) 



which both agree in the Landau gauge and the variables in each expression correspond to those 
of the scheme indicated in the left hand side. 



4 Tensor current in the RI' scheme. 

We now turn to the computation of the anomalous dimension for the flavour non-singlet tensor 
current at zero momentum insertion in the Landau gauge. This calculation is similar to the 
one for the vector current though the decomposition of the Green's function G'^^iu{p) will have 
different Lorentz structures. In particular we have 



(4.1) 

since cr'''^ is antisymmetric in its Lorentz indices. The components are deduced from 



.(2) 



.(1) 



ip) 



1 



.(2) 



•^Ap) = - 



4(d-l)(d-2) 
1 

4(d-l)(d-2) 



tr 



(tTM-Gj^^.^(p)) + ^2^^ [i^j^Pu - ^luP,.)G^"^,.^ip)) 
tr (a^^Gj;^,^(p)) + Atr((^7^p,-|^7,p^)Gj;^.^(p)) 



(4.2) 



20 



As the anomalous dimension of the current had been computed for arbitrary covariant gauge 
parameter in the MS scheme, |351 1361 1^ . we note that the definition of the RI' scheme renor- 
mahzation constant is 



hm 



(1) 



1 . 



(4.3) 



Therefore, we find that the gauge dependent renormahzation constant is 
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where we have used the same symbohc manipulation programme to compute this as for the 
scalar and vector cases aside from changing the Feynman rule for the operator insertion. More- 
over, given that the programme correctly reproduces the gauge independent MS renormalization 
constant for all a, [351 1361 151] . we are confident that H4.4|) is correct. Therefore, from 
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we find that, in four dimensions, 
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where the four loop expression of ()4.7|) is available for QED in the quenched approximation, 
pS] . Specifying the Landau gauge for the colour group 5C/(3) we find 
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1/2, Cp = 4/3 and Ca = 3. Finally, we note 
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which implies 
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For S^^) and S^^^ 



it transpires that there are no contri- 



butions to the finite part to and including three loops. This is consistent with the one loop 
evaluation of the same quantity in |37j when the chiral limit is taken. For completeness we have 
again checked our RI' scheme anomalous dimension by constructing the conversion function 
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we again have exact agreement. 



5 Discussion. 



To conclude we have first renormahzed QCD to three loops in arbitrary covariant gauge in 
the RI' scheme. The full renormalization was necessary since, for example, the anomalous 
dimension of the gauge parameter is required when converting the renormalization constants 
to renormalization group functions for non-zero a. Although in practice one only requires 
information in the Landau gauge, computing for a 7^ provides important internal checks on 
the calculation such as for comparing with established MS results. Further, we have extended 
the machinery of ^2 to compute the anomalous dimensions of the tensor operator ipcT^^^ip at 
three loops in the chiral limit in the Landau gauge in a scheme which is natural in lattice 
regularization. Given this approach it would be interesting to examine other operators whose 
anomalous dimensions are required in the RI' scheme such as the low moments of the twist-2 
Wilson operators which occur in the operator product expansion in deep inelastic scattering and 
those relating to transversity, in order to provide the foundation to improve lattice estimates of 
matrix elements. 
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